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Abstract 

A coarse-grained model for solutions of polymers in supercritical flu- 
ids is introduced and applied to the system of hexadecane and carbon 
dioxide as a representative example. Fitting parameters of the model to 
the gas-liquid critical point properties of the pure systems, and allow- 
ing for a suitably chosen parameter that describes the deviation from the 
Lorentz-Berthelot mixing rule, we model the liquid-gas and fluid-fluid un- 
mixing transitions of this system over a wide range of temperatures and 
pressures in reasonable agreement with experiment. Interfaces between 
the polymer-rich phase and the gas can be studied both at temperatures 
above and below the end point of the triple line where liquid and vapor 
carbon dioxide and the polymer rich phase coexist. In the first case in- 
terfacial adsorption of fluid carbon dioxide can be demonstrated. Our 
model can also be used to simulate quenches from the one-phase to the 
two-phase region. A short animation and a series of snapshots help to 
visualize the early stages of bubble nucleation and spinodal decomposi- 
tion. Furthermore we discuss deviations from classical nucleation theory 
for small nuclei. 
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1 Introduction 



The phase behavior of polymer-solvent systems has important application in the 
industry for the production and processing of many kinds of plastic materials ^ 
I2»i3,^^£|. An example is the formation of solid polystyrene foams EI El El- 
In addition, both the understanding of the equilibrium phase diagram of these 
systems and the kinetic mechanism of phase separation are challenging problems 
of statistical mechanics. On the one hand, variation of the molecular weight of 
the polymer (i.e. the "chain length" of the flexible linear macromolecule, N) 
offers a control parameter that leaves intermolecular forces invariant, and hence 
allows a more stringent test of the theory than would be possible for small 
molecule systems. On the other hand, it is crucial to take the compressibility 
of the system fully into account, since for practical reasons one normally wishes 
to use supercritical fluids (e.g., CO2) as a solvent (Ij i2j ^ ,4, ,5j i6j . Then small 
changes of the pressure result in a large variation of the solvent density, and 
this fact obviously facilitates applications. While for an incompressible polymer 
solution (as it is modelled, for instance, by the well-known Flory-Huggins lattice 
theory (TIEIEIEI) the early stages of phase separation can only be studied if one 
realizes a rapid temperature quench from a state in the one phase region into the 
miscibility gap (cf. Fig. for compressible polymer solutions it is possible to 
induce phase separation by an (experimentally much more convenient) pressure 
jump (cf. Fig. n^)- Of course, in Fig. ^ we have made the tacit assumption 
that including the pressure does not change the qualitative form of the phase 
diagram, i.e., the point {Tc,Xc) in (T, x,p)-space becomes a line {Tc{p),Xc(j})) 
or, expressed in the appropriate inverse functions, (j)c{T), Xc{T)), respectively 
(Fig. ^p). In reality, the situation is more complicated, since one must treat 
both the densities of CO2 and hexadecane as two coupled order parameters, 
which both play a crucial role in the possible phase transitions (remember that 
for the pure solvent at its critical point gas- liquid phase separation starts). As a 
result, it can happen that both gas-liquid and liquid- liquid phase separation are 
present in the phase diagram, leading to the occurrence of a triple line (Fig. [SJ 
As will be discussed below, the presence of the triple line has a profound 
influence on the interface properties and hence the nucleation behavior in the 
system. 

The present paper is devoted to a study of the early stages of phase separa- 
tion kinetics in such systems by computer simulation. We use a coarse-grained 
model for mixtures of hexadecane and carbondioxide as an archetypical exam- 
ple. The reason for the choice of this particular system is the fact that recently 
extensive experiments on nucleation in this system have been performed 12 . In 
Sec. 2, we shall briefly describe the model and the simulation techniques. Sec. 3 
presents the relevant static equilibrium properties, in particular, the phase dia- 
gram and the properties of planar interfaces separating coexisting phases. Sec. 4 
describes simulations of quenching experiments under different conditions, that 
lead either to nucleation and growth or spinodal decomposition. Sec. 5 summa- 
rizes our conclusions. 
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2 Model and Simulation Techniques 



A fully atomistic simulation of mixtures of hexadecane (Cig H34) and carbon- 
dioxide (CO2), which aims at both establishing phase diagrams as a function of 
the three variables T,p and x, and the study of phase transition kinetics under 
various conditions, would be an extremely formidable problem. The potentials 
for the length of the covalent bonds in these molecules (as well as the potentials 
for bond angles) are very stiff. Hence, an extremely short time step (of the 
order of about 1 fs) would be required in a Molecular Dynamics (MD) simula- 
tion. In a corresponding Monte Carlo (MC) simulation only very small random 
displacements of atoms would be admissible. Thus, for liquid alkanes it is a 
quite common and well-established practice to integrate CH2-groups (as well as 
the CH3 end groups) into "united atoms" and also to work with a fixed 
C-C bond length. These simplifications already reduce the necessary computer 
time by a factor of 100. However, we have estimated that even for such a sim- 
plified model computer resources as they are available today still are not yet 
sufficient. Thus it was chosen to simplify the model even further, representing 
CO2 by a single pseudo-atom, and representing Cig H34 by a flexible chain of 5 
subsequent effective segments, each of which then contains roughly 3 successive 
C-C bonds. All such effective monomeric units interact with a truncated and 
shifted Lennard Jones (LJ) potential, 

VLj{r) = 4e,.^A{ — r - { — )' + T^], r < = 2 . 2'/'a , (1) 
r r 1d384 

while VLj(r > Tc) — 0. Subsequent effective monomers along a chain are in 
addition exposed to a finitely extensible nonlinear elastic (FENE) potential 

I^FENE(r) = -33.75ew,ln[l - (7^^)'] (2) 

The constants in this potential are chosen such that the most favorable distance 
between bonded neighbors is 0.96 CThh, while the preferred distance between non- 
bonded effective monomers is 2^/^ Uhh ~ 1-12 CThh- This mismatch is desirable 
to prevent crystallization at high densities which is appropriate for glass- 
forming polymers |31 E El E] • 

The parameters CThh, Chh set the scales for energy and length of our hexade- 
cane model. Since we wish that our model fits the thermodynamic properties 
of this material in the fluid phase as faithfully as possible, we have adjusted 
them such that the critical temperature — 723 K and the critical density 
Pc — 0.21 g/cm^ of hexadecane are correctly reproduced. In simulations, a 
critical point can be computed by a finite size scaling study of the liquid-gas 
transition along the Hues of the techniques proposed by Wilding et al. ^Hl- A 
comparison of critical temperature and density in Lennard- Jones units with the 
experimental values yields Ehh = 5.787 • lO'^i J, CThh = 4.523 • 10-i°m [TTIIT^ . 

Since the description of hexadecane is thus already reduced to a crudely 
coarse-grained model, it would not make sense to keep all atomistic detail for 
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carbon dioxide. Thus the CO2 molecule is also coarse-grained into a point par- 
ticle, and we require an interaction potential between CO2 molecules of exactly 
the same LJ form as in Eq. J^l, but with parameters £cc, <^cc- Requiring once 
more that the critical temperature and density of real CO2 are correctly re- 
produced we obtained [I7| £cc = 4.201 • IQ-^^ J, and CTcc = 3.693 • 10-^° m. 
It is clear that this procedure ignores some physical effects, e.g., the fact that 
the CO2 molecules carry electrical multipole moments is completely neglected. 
Nevertheless our model not only describes the gas-liquid coexistence curve of 
hexadecane and CO2 over a reasonable temperature range J7|, but also other 
data (like the critical pressure, or the temperature dependence of the surface 
tension near Tc |18p. 

Of course, the choice of interaction parameters between CO2 and hexadecane 
is more subtle. We model the interaction between the pseudoatoms representing 
CO2 and the pseudoatoms representing three subsequent CH2 groups again by 
a LJ potential of the form Eq. flj, but now with parameters £hc, o-\\c- It then 
remains to find an optimal choice for these parameters. The simplest choice 
would be the well-known Lorentz-Berthelot mixing rule |19j 

£hc ~ \J £hh Ecc , Che (Chh + 0-cc)/2 ■ (3) 

However, when one tries the choice Eq. Q the resulting phase diagram of our 
model system would be of type I in the classification scheme of Van Konynen- 
burg and Scott [20112] (i-e., the critical points of pure hexadecane and CO2 are 
connected by a critical line of the mixture system. Liquid-liquid phase seper- 
ation does not exist and thus there is no three phase coexistence). However, 
experimentally it is known [23 that alkane - CO2 mixtures exhibit a type I - 
phase behavior only for very short alkanes, while for hexadecane - CO2 mix- 
tures the phase diagram is of type HL Instead of a connecting line we rather 
observe a topology as shown in Fig. |21 El ESI- In the (p,T) projection of the 
phase diagram, a critical line emerges from the critical point of hexadecane and 
does not end at the critical point of CO2. In addition we observe liquid-liquid 
immisciblity and a three-phase line. When we want to obtain this behavior from 
our coarse-grained model we must allow for a deviation from the first equation 
of Eq. (O , assuming rather 

Ehc = ^V^hh, £cc , ■ (4) 

It turns out that a suitable choice for the parameter ^ which characterizes the 
deviation from the Lorentz-Berthelot mixing rule is |17[ 1181 123 C = 0.886. 
Eq. Q with this choice of ^ was in fact used to compute the phase diagram 
shown in Fig.|31 However, experimental measurements of the critical line vary 
considerably and a small modification in ^ rises or lowers the critical line. Details 
about how such phase diagrams are in fact estimated from the simulation are 
given elsewhere |18| . 

For a study of the phase diagram (Fig. (JJ as well as for a study of the 
interfacial free energy between coexisting gas and liquid phases ^1 E| the 
grandcanonical ensemble is used, where volume V ^ temperature T, and the 
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chemical potentials fj,c, fJ-h are fixed. Both the particle numbers Nc and iVh 
of CO2 and hexadecane and the pressure p are then (fluctuating) observables 
"measured" in the simulation ^]. Fig.2|shows an isothermal slice through the 
phase diagram at T = 486 K. This temperature is higher than the temperature 
of the critical end point where the critical line ending at the critical point of 
CO2 and the line of triple points {ptiip{T)) in the (p, T) phase diagram (Fig.|3Jl 
meet. Therefore unlike Fig.|2]there is no triple point in the phase diagram, and 
the qualitative features are the same as in Fig.QJj. At this temperature we shall 
examine the kinetics of the phase separation in Sec. 4. 

In order to observe the kinetics of bubble nucleation in the vicinity of the 
binodal we control the undersaturation of the system by fixing the chemical po- 
tential of both species. As a starting configuration we use a homogeneous state 
equilibrated at a higher temperature. This situation corresponds to a simulation 
cell in contact with a much larger reservoir which is held at constant undersat- 
uration. As the number of particles in the simulation cell is not conserved, and 
the particles do not move according to a realistic dynamics, we do not obtain 
information about the time scale of bubble nucleation. The kinetics of phase 
separation in the vicinity of the binodal is, however, chiefly determined by the 
free energy barrier the system encounters on its path towards the equilibrium 
state. Hence, we expect the relaxation path to be similar to a simulation with 
realistic dynamics. The most natural choice would be to apply Molecular Dy- 
namics methods |^|^. However, if the number of particles was conserved we 
would need prohibitively large simulation cells for the undersaturation not to 
decrease substantially even in the very early stages of nucleation. The simula- 
tion of spinodal decomposition took place in the canonical ensemble. Again, a 
homogeneous starting configuration was equilibrated at a higher temperature 
and quenched to a state below a critical point of the mixture system. We only 
allowed for local Monte Carlo displacements which should yield kinetics compa- 
rable to those obtained in Molecular Dynamics. 

3 Static equilibrium properties 

First we return to the phase diagrams in Figs. and discuss their properties 
more closely. Note that at the temperature chosen in Fig. 0] (T = 486 K) we 
are very far below the critical point of pure hexadecane (which is = 723 K). 
Therefore the hexadecane melt (with some dissolved carbon dioxide) coexists 
with almost pure CO2 the composition of which gradually changes from gas- to 
liquid-like with increasing pressure. In both cases there is almost no hexadecane 
present. This is reflected in the coexistence curve which rapidly approaches the 
CO2 molar fraction x ^ 1. 

We have also included the results of an analytical approximation for the 
phase diagram of our model, obtained from a thermodynamic perturbation the- 
ory [231 called TPTl, developed along the lines of Wertheim et al. [27|. It is seen 
that this theory describes the coexistence curve in Fig0]at pressures that are 
much smaller than the critical pressure Pcrit {T) almost quantitatively, while the 
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critical pressure itself is clearly overestiraated. Of course, such an overestima- 
tion of the critical parameters Pcrit (T) , Xcrit (T) is quite typical for all theories 
that invoke a mean field description of the critical behavior, as TPTl does. A 
similar analytical study for more alcane-C02 systems has been performed re- 
cently I2H1- Note also that the TPTl theory readily yields a spinodal curve, 
which has the standard meaning of separating "metastable states" from "unsta- 
ble states" in the phase diagram |29II30[I^ . The naive (mean- field) description 
implies that in the metastable region of the phase diagram the initial stages of 
phase separation kinetics are described by nucleation and growth [221 EOl ; while 
in the unstable region the decay mechanism is spinodal decomposition |3U[ 13 1| . 
In systems with short range forces there is no well-defined sharp spinodal line, 
the actual transition between both decay mechanisms occurs rather gradually 
in a broad transition region, and this region is not centered around the mean 
field spindodal but occurs close to the coexistence curve However, for 

the quenching experiment performed at T = 486 K and x = 0.60, where the 
final state corresponds to a pressure of about p «130 bar, we are so close to 
the coexistence curve and so far from the mean field spinodal, that classical 
nucleation- type behavior |29[l30j should be observable. In order to demonstrate 
the mechanism of spinodal decomposition, we quenched into the unstable region 
as defined by the mean- field spinodal (compare with Fig. 0)). 

Classical nucleation theory assumes (in our case) the formation of spher- 
ical gas bubbles of essentially pure CO2 within a hexadecane matrix . For 
a quantitative understanding of the free energy barrier against homogeneous 
nucleation, clearly the understanding of the interface between macroscopic co- 
existing phases is a prerequisite. As discussed elsewhere ^|^], such interfaces 
are most conveniently studied in the grandcanonical ensemble in conjunction 
with a multicanonical preweighting scheme that generates mixed-phase config- 
urations. Choosing a simulation box of rectangular shape Ly.Ly.cL with 
a constant c = 2 or larger and periodic boundary conditions, one generates 
system configurations where the two coexisting phases are separated by two 
parallel L y L interfaces (compare with Fig. ^C). "Measuring" the canonical 
probability of these two-phase configurations relative to the probability of the 
pure phases is a standard method for the estimation of interfacial free energies 
[TB1[T7,"TF,15",1!F,'M',^,36 , 37 . The results for T=486 K are shown in the in- 
let of Fig.^ In addition, this method can be used to generate well-equilibrated 
configurations of interfaces allowing a study of the interfacial structure. Fig. El 
shows, as an example, two snapshot pictures of such states containing inter- 
faces, one at T=486 K and the other at T=243 K. The size of the particles is 
enlarged (radii shown are CTcc, CThh rather than crcc/2, (Thh/2), for the sake of 
a clearer view. The temperature T=243 K is slightly lower than the tempera- 
ture of the critical point of pure CO2, and thus the system possesses a triple 
point where a polymer-rich phase coexists with two C02-rich phases, one being 
liquid CO2 and the other the gas. Consequently, when we study the interface 
between the polymer-rich phase and gaseous CO2, a layer of the third phase 
(liquid CO2) intrudes at this interface. The thickness of this layer is expected 
to diverge when the pressure approaches the value of the triple point pressure. 
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For higher pressures, we then have a coexistence between the polymer rich- 
phase and the hquid-Hke dense CO2, similar to the situation that occurs at the 
higher temperature (Fig.|SJi). At low enough temperatures, where a triple point 
occurs, one expects a strong decrease of the interfacial tension already when 
the triple point pressure is approached : the interfacial tension between the 
polymer-rich phase and liquid-like CO2 is much smaller than the interfacial 
tension between the polymer-rich phase and the gas. The barrier against nu- 
cleation also decreases when one approaches the spinodal from the metastable 
region. This fact facilitates nucleation of relatively large droplets observable on 
the time-scales of our Monte Carlo simulation. 

4 Monte Carlo simulation of bubble nucleation 

The animation in Fig. |^ visualizes how phase separation in polymer solutions 
proceeds via nucleation. In the beginning the system fluctuates around a meta- 
stable free energy minimum. One observes the "birth and death" of small den- 
sity fluctuations in the hexadecane matrix (displayed as red spheres) which 
become visible whenever the blue background shines through the slice. These 
irregular voids are too small, however, to lead to an immediate decay of the 
initial metastable state. Usually the voids also contain a few CO2 molecules 
(displayed as small yellow spheres). Only after some time lag a void manages 
to grow to critical size and overcome the free energy barrier which seperates the 
metastable from the homogeneous equilibrium state. From now on the bubble 
will grow until it fills the whole simulation box. As expected, the critical void 
gets also filled with CO2 molecules, thus decreasing the supersaturation of the 
remaining polymer-rich phase. Remarkably, this filling of the bubble does not 
occur homogeneously: at the interface between the gas and the polymer, i.e., 
at the surface of the bubble, the CO2 density is clearly higher than it is in the 
center of the bubble. This surface enrichment effect may be interpreted as a 
precursor to interfacial wetting of the fluid CO2 predicted when one approaches 
the triple point . A quantitative analysis of this phenomena is in preparation 

SOI. 

Of course, the observation of a nucleation event displayed in Fig. is by no 
means the first observation of nucleation phenomena by computer simulations: 
e.g., already Ref. [2^1 contains a series of snapshot pictures illustrating nucle- 
ation in the two-dimensional lattice gas model. Since then many more elaborate 
studies have appeared. However, the present work deals with nucleation in a 
system of industrial relevance, namely bubble nucleation in metastable poly- 
mer melts supersaturated with supercritical CO2. In addition, the nucleation 
mechanism here is quite nontrivial, since the "critical bubble" (which just has 
the size to be at the nucleation barrier, where there is a 50% chance of growth 
or decay) is characterized by two variables, its size and the number of CO2 
molecules it contains. Similarly the interface is described by nontrivial profiles 
of two variables, the total density and the relative concentration of CO2. Such 
a nontrivial coupling between two order parameters has also been predicted by 
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recent self-consistent field calculations based on a qualitatively similar model 
Since these mean field descriptions do not take into account the effects of 
fluctuations, a check by computer simulation methods is clearly warranted. 

When we quench to a position below the critical point (Fig. [7^, the be- 
havior is rather different: fluctuations in the relative concentration of CO2 get 
gradually more pronounced everywhere. These fluctuations are not localized, 
however, in the form of identifiable bubbles, but rather form an irregular perco- 
lating network. This is the hallmark of spinodal decomposition |3(JI 1311 . Since 
the linear theory of spinodal decomposition only holds for systems with long 
range forces in the very early stages we do not attempt a more quantita- 
tive analysis of this data. Note also that the build-up of these concentration 
fluctuations is relatively faster than the nucleation of bubbles. Hence the sep- 
aration of time scales between the structural relaxation and the rate constant 
of spinodal decomposition is less well established. In the intermediate stages 
of spinodal decomposition, we also expect hydrodynamic mechanisms to prevail 
Pnil^, which are not captured by Monte Carlo moves. 

Returning now to the simulation of nucleation phenomena, we ask how one 
can quantify these observations. A straightforward but tedious method would 
be to repeat dynamic simulations as shown in Fig. |^ many times in order to 
extract an estimate of the nucleation rate from the average time it takes to 
form a bubble or a droplet. Of course, since the nucleation rate is expected to 
decrease dramatically if the quench depth is decreased slightly, such a direct 
method works hardly in practice. An indirect method to obtain information 
on the surface free energy of clusters as a function of cluster size simply uses 
the final equilibrium states of simulations. An estimate for the free energy as a 
function of density can be obtained directly from the probability distribution. In 
this context special care has to be taken of finite-size effects to identify regions of 
densities which correspond to a single cluster |41[I42| . If we assume that bubbles 
(or droplets) are spherical, filled with gas and surrounded by a homogeneous 
liquid of coexistence density (or vice versa), we can assign a radius to each 
density and hence a free energy estimate to each radius. Fig. |S1 indicates that 
the surface free energy obtained this way is somewhat smaller than expected 
from the so-called "capillarity approximation" . A simple estimate of the free 
energy is Fg = 47ri?^7int, where R is the radius of the cluster and 7int the 
interfacial tension of flat interfaces as shown in Fig. [3 and derived in Fig. 01 
As a consequence, the nucleation barrier is lower than the classical theory of 
homogeneous nucleation [231201 would predict. Differences in 7int also decrease 
with increasing cluster size and distance to the critical point of the mixture. At 
this point, however, the findings are still somewhat preliminary, and more work 
will be necessary to obtain sufficiently accurate results |43j . 

5 Concluding Remarks 

In summary, we have presented a simple model and computational tools which 
allow us to examine equilibrium properties and kinetics of phase separation 
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in systems of industrial relevance. We were able to produce a two-component 
mixture phase diagram with a complete critical line which is in qualitative agree- 
ment with experiments. For a system at T=486 K we determined the isotherm 
and interfacial tension of a flat surface as a function of pressure. In a small ani- 
mations and a series of snapshots we visualized bubble nucleation and spinodal 
decomposition for this specific system. For the first case we found a precursor 
of interfacial wetting which was also observed on a flat surface in the vicinity 
of the three-phase coexistence line. Finally, we found evidence that for small 
bubbles the interfacial tension is smaller than expected by classical theory of 
homogeneous nucleation. 
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Figure 1: a) Schematic phase diagram of an incompressible polymer solution. 
Temperature T and molar fraction x of the solvent are variables, pressure p is 
constant. The position of the critical point is shown by •. We indicate how a 
sudden decrease of temperature from a state in the one phase region to a state in 
the two-phase region is carried out. Note that the molar fraction of the critical 
point tends to unity when the chain length N of the macromolecule gets large, 
b) Schematic isothermal slice of a phase diagram of a compressible polymer 
solution, using pressure p and molar fraction x of the solvent as variables. We 
indicate how a pressure jump experiment could be performed. In both cases 
a) and b) it is assumed that the state after the jump is in the unstable region 
of the phase diagram, i.e., underneath the spinodal curve. If the quench ends 
in between the spinodal curve and the coexistence curve, phase separation will 
start by homogeneous nucleation. 
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Figure 2: Isotherms of a simple model-system obtained from self-consistent 
field calculations 11.. At T=343.7 K the phase diagram is like the schematic 
diagram in Fig. ^p. However, at temperatures close to the critical point of pure 
CO2 we observe liquid-liquid immiscibility for a type III system. A three phase 
line and a second critical point occur. 
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Figure 3: Projection of the hexadecane-C02 phase diagram (p,T,x) onto the 
pressure-temperature plane. Black lines, + and •: simulations (for ^ — 0.886), 
red lines: analytic calculations from a thermodynamic pertubation theory 
(TPTl) based on the same model [23: blue lines: experimental data. Gas- liquid 
coexistence lines are from 0S] , the two critical lines are taken from 39: and 
| 22l fIH| (dashed blue) and differ considerably. Experiments are in qualitative 
agreement with both simulations and analytic calculations. Three features can 
be identified: Gas-liquid coexistence lines {+) of pure CO2 and hexadecane both 
end at the corresponding critical point (at T=304 K and T=723 K respectively). 
A line of critical points (•) emerges from the critical point of pure hexadecane 
and gradually changes its composition from gas-liquid hexadecane to liquid CO2 
- liquid hexadecane. The red dotted three-phase line (TPTl) lies slightly below 
the corresponding CO2 coexistence curve, but only above the critical point of 
CO2 it is distinguished in this plot. The critical point of CO2 and the the end 
point of the three-phase line are connected by another short critical line which 
cannot be distinguished, too. 




Figure 4: Isothermal slice through the phase diagram in Fig. U] at temperature 
T=486 K. Black curve with circles shows the simulated coexistence curve. Red 
curves are corresponding predictions from TPTl theory for coexistence 
(full) and the spinodal curve (dotted) , respectively. The position of the critical 
points are denoted by •. o indicate phase diagram positions of nucleation and 
spinodal decomposition movies in Fig. Eland Fig.[7| (Inset): surface tension as 
a function of pressure. 
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Figure 5: Snapshot pictures of systems in the center of the two-phase coexis- 
tence region for T=243 K and T=486 K. Snapshots show a 2ahh slice through 
a box with dimensions ^ Ly x =18 Ohh x 18 Ohh x 54 ahh- Positions 
of particles are projected into the xz plane. CO2 particles are shown as yellow 
circles of radius CTcc, while effective monomers of hexadecane are shown as red 
circles of radius (Thh (we choose dhh = 1 as our unit of length, which then implies 
(Tec = 0.816), the background is blue. For T=243 K interfacial wetting can be 
observed. 



Figure 6: (Remark: The final version (submitted to New Journal of Physics) 
will contain a short animation. Here, only snap-shot pictures are shown.) A 
short animation showing the time evolution of a quenching experiment into the 
metastable region (r=486 K, x=0.60 and (final) pressure p w 130 bar - compare 
with FigQ. Phase separation occurs via nucleation. The linear dimension of the 
simulation box is 22.5 ahu- The movie is generated from subsequent snapshot 
pictures of a grandcanonical Monte Carlo simulation. For clarity, not the whole 
simulation but only a slice of thickness 2 auh is shown. The color coding is the 
same as in Fig. [S] 



Figure 7: A series of snapshots visualizing the early stages of spinodal de- 
composition after a quenching experiment into the unstable region (T=486 K, 
x=0.833 and (final) pressure p « 52 bar - compare with Fig. Phase separa- 
tion occurs via spinodal decomposition. The size of the box is 91.21 (Jhh x 91.21 
ahh X 17.53 ahh, but only a slice of 10 <7hh is shown. The system is actually 
more dilute than the one shown in Fig. |B| (a) Starting configuration, (b) after 
125000 and (c) after 500000 local displacements (0-0.3 0-5) per particle. The 
color coding is the same as in Fig. |31 
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Figure 8: Free energy as a function of droplet size for T=486 K and p = 0.5 Pcrit 
and pwO (coexistence pressure of pure hexadecane - compare with Fig. |31 and 
Fig. EJ. Indigo line: A simple estimate for the free energy: A_F — An'yR'^, 
surface tension 7 is taken from Fig. 0] (inlet) (flat surface). Black: results 
from simulations of different system sizes. Only the envelope of the curves is 
relevant. Other parts of the curves belong to regions of the distribution where 
no nucleation is expected. For small droplets the free energy is smaller than 
expected. Differences decrease with increasing radius and decreasing distance 
from the critical point of the isotherm. 
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